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ABSTRACT 

We report the detection of spatially extended CO( J=l— >0) and CO( J=5— >A) emission in the z=2.49 
submillimeter galaxy (SMG) J123707+6214, using the Expanded Very Large Array and the Plateau 
de Bure Interferometer. The large molecular gas reservoir is spatially resolved into two CO( J=l— >0) 
components (north-east and south-west; previously identified in CO J=3— >2 emission) with gas masses 
of 4.3 and 3.5 xlO 10 (aco/0-8)M Q . We thus find that the optically invisible north-east component 
slightly dominates the gas mass in this system. The total molecular gas mass derived from the 
CO(J=l— >0) observations is >2.5x larger than estimated from CO(J=3— >2). The two components 
are at approximately the same redshift, but separated by ^20kpc in projection. The morphology 
is consistent with that of an early-stage merger. The total amount of molecular gas is sufficient to 
maintain the intense 500M©yr _1 starburst in this system for at least ^160 Myr. We derive line 
brightness temperature ratios of r 31 =0.39±0.09 and 0.37±0.10, and r 5 i=0.26±0.07 and 0.25±0.08 
in the two components, respectively, suggesting that the J>3 lines are substantially subthermally 
excited. This also suggests comparable conditions for star formation in both components. Given the 
similar gas masses of both components, this is consistent with the comparable starburst strengths 
observed in the radio continuum emission. Our findings are consistent with other recent studies that 
find evidence for lower CO excitation in SMGs than in high-z quasar host galaxies with comparable 
gas masses. This may provide supporting evidence that both populations correspond to different 
evolutionary stages in the formation of massive galaxies. 

Subject headings: galaxies: active — galaxies: starburst — galaxies: formation — galaxies: high- 
redshift — cosmology: observations — radio lines: galaxies 



1. INTRODUCTION 

Detailed studies o f subm illimeter galaxies (SMGs; see 
review by Blain et al..20Q2j) have revealed that they repre- 
sent a relatively rare, but cosmologically important high 
redshift population of massive galaxies. They harbor in- 
tense (>500 Mq yr _1 ), often heavily obscured, but rather 
short-lived (<100Myr) starbursts that rapidly consume 
their gas content through star formation at high efficien- 
cies. SMGs may trace a common phase in the formation 
and evolution of massive galaxies in the early universe, 
making them the likely progenitors of today's massive 
spheroidal galaxies. 

Given their substantial dust obscuration, the most in- 
sightful way to study SMGs and their star formation 
properties is through the dust-reprocessed emission at 
rest-frame far-infrared (FIR) wavelengths. While contin- 
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uum diagnostics at such wavelengths are particularly use- 
ful to determine star formation rates (SFRs) unaffected 
by obscuration, the most insightful way to study the fate 
of such galaxies is through emission line fluxes, morphol- 
ogy and dynamics of the material that fuels the star for- 
mation, i.e., molecular gas (typically CO). Molecular gas 
was detected in >30 SMGs to date, revealing large gas 
reservoirs of >10 1 M Q in most cases (see Solomon & 
Vanden Bout [2005| for a review). However, most of these 
studies were carried out in mid- to high- J CO transitions, 
rather than the fundamental CO(J=l— >0) transition. 

Recent studies of CO(J=l— >0) emission indicate that 
(in contrast to high-z quasar hosts) this line appears to 
carry a higher brightness temperature than the mid- to 
high-J lines in several SMGs, sug gestin g relatively low 
gas excitat ion (e.g., Hainline et al. 120061 : R iecher s et al. 
120061 120101: Carilli et al. 120101 : Ivison et al. 120101 Harris 
et al. 2010). Therefore, earlier studies may systemati- 
cally underestimate the total amount of molecular gas 
that is present in SMGs. The low excitation in at least 
part of the molecular gas reservoirs raises the question 
whether the commonly used aco conversion factor from 
CO luminosity to gas mass (Af gas ) for ultra-luminous in- 
frared galax ies (U LIRGs) in the nearby universe (Downes 
& Solomon I1998D is applicable, or if this practice leads 
to an underprediction of Af gas - 

To further investigate this issue, we have started a 
systematic Expanded Very Large Array (EVLA) survey 
of CO(J=l— >0) emission in high-z SMGs and quasar 
host galaxies that were previously detected in higher- 
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Fig. 1.— EVLA CO(J=1->0) emission toward J123707+6214. Left: Velocity-integrated CO(J=1->0) map over 490 kms -1 (54MHz). 
At a resolution of 2.8"x2.1" (as indicated in the bottom left), the emission is resolved into a n orth-e ast (ne) and south-west (sw) com- 
ponent. The crosses indicate the peak positions of the CO(J=3— >2) emission (Tacconi et al. 120061) . Contours are shown in steps of 
lcr=39 fi]y beam - 1 , starting at ±3o\ Right: Beam-corrected CO(J=l— >0) spectra toward the ne (top) and sw (bottom) components at 
18 kms" 1 (2 MHz) resolution (histograms), along with Gaussian fits to the line emission (black curves). The velocity scale is relative to 
the tuning frequency of 33.029 GHz (corresponding to ,2=2.49). 



J CO lines. In this Letter, we report the first re- 
sults from this study, i.e., the detection of spatially ex- 
tended CO(J=1^0) emission toward the z =2.49 SMG 
J123707+6214 (GN19; HDF242; Borys et al. 120031) . We 
also report the detection of CO(J=5— >A) emission in 
this system, using the Plateau de Bure Interferometer 
(PdBI). We use a concordance, flat ACDM cosmology 
throughout, with g = 71 kms^M pc" 1 , ft M =0.27, and 
fi A =0.73 (Spergel et al. 1200311200^ . 



of overlapping edge channels between sub-bands; corre- 
sponding to ~2300 kms -1 at 9.08mm) with the WIDAR 
correlator. 

For data reduction and analysis, the AIPS package was 
used. All data were mapped using 'natural' weighting. 
Maps of the velocity-integrated CO J=1^0 line emis- 
sion yield a synthesized clean beam size of 2.8" x 2.1" 
at an rms noise level of 39 /iJy beam -1 over 490 kms -1 
(54 MHz). 



2. OBSERVATIONS 
2.1. EVLA 

We observed the CO(J=1^0) (u Iest = 115.2712GHz) 
emission line toward J123707+6214 using the EVLA. At 
z=2.49, this line is redshifted to 33.029 GHz (9.08 mm). 
Observations were carried out under excellent weather 
conditions (typical atmospheric phase rms: 2.3° on a 
300m baseline) in D array on 2010 April 11 (NRAO 
Legacy ID: AR708), resulting in 2.9 hr on-source time 
with 18 antennas (equivalent to 1.2 hr with 27 antennas) 
after rejection of bad data. The nearby (5.4° distance) 
quasar J1302+5748 was observed every 7.5 minutes for 
pointing, secondary amplitude and phase calibration. 
For primary flux calibration, the standard calibrator 
3C286 was observed, leading to a calibration that is accu- 
rate within <10%. Observations were set up using a total 
bandwidth of 252 MHz (dual polarization; after rejection 



2.2. PdBI 

We observed the CO(J=5-^4) (i/ rest = 576.2679 GHz) 
emission line toward J123707+6214 using the IRAM 
PdBI. At z=2.49, this line is redshifted to 165.1198 GHz 
(1.82mm). Observations were carried out under good 
weather conditions with five antennas in the compact 
D configuration during 3 tracks on 2008 July 10, and 
August 09 and 14 (IRAM program ID: SC47), for a 
total of 14.9 hr, resulting in 8.1 hr of 6 antenna equiv- 
alent on-source time after rejection of bad data. The 
nearby quasars B0954+658 and B1418+546 (distance 
to J123707+6214: 17.4° and 15.4°) were observed ev- 
ery 20 minutes for pointing, secondary amplitude and 
phase calibration. For primary flux calibration, the stan- 
dard calibrators MWC 349 and 3C454.3 were observed. 
Observations were set up using a total spectrometer 
bandwidth of 1 GHz (dual polarization; corresponding 
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Fig. 2.— PdBI CO(J=5->4) emission toward J123707+6214. Left/Middle: Velocity-integrated CO(J=5-s>4) maps over 563/309 kms -1 
(320/180 MHz). At a resolution of 3.7"x3.0" (as indicated in the bottom left), the emission is resolved into the ne and sw components 
also seen in CO(J=l— >0). The velocity ranges are selected to maximize the integrated (left) and sw component (middle) signal-to-noise 
ratio. The crosses are the same as in Fig.[T] Contours are shown in steps of 1<t=0. 39/0. 54 mjy beam -1 , starting at ±2o\ Right: Integrated 
CO(J=5— >4) spectrum at 36 kms -1 (20 MHz) resolution (histogram), along with a Gaussian fit to the line emission (black curve). The 
velocity ranges of the two maps are indicated by the black lines. The velocity scale is relative to the tuning frequency of 165.1198 GHz 
(corresponding to £=2.49). 



to -1800 kms -1 at 1.82mm). 

For data reduction and analysis, the GILDAS pack- 
age was used. All data were mapped using 'natural' 
weighting. Maps of the velocity-integrated CO J=5— s4 
line emission yield a synthesized clean beam size of 
3.7" x 3.0" at an rms noise level of 0.39/0.54 mjy beam -1 



over 563/309 kms 



(320/180 MHz). 
3. RESULTS 



3.1. Gas Morphology and Emission Line Properties 

We have detected spatially resolved CO(J=l— »0) line 
emission toward the z=2.49 SMG J123707+6214, mea- 
suring two components at la (north-east; 'ne'jfl and 6a 
(south-west; 'sw') significance in the velocity-integrated 
emission line map (Fig. [IJ left). We do not detect the 
underlying continuum emission at a 3<r upper limit of 
54/iJybeam -1 at 9.08mm (rest-frame 2.6mm). 

From Gaussian fitting to the CO( J=l— s>0) line profiles 
of the ne and sw components (Fig.[lJ right) F^l we obtain 
line peak strengths of S'„=375±57 and 340±64 /iJy at line 
FWHMs of du=454±87 and 414±92kms -1 , centered 
at redshifts of z=2.4879±0.0005 and 2.4873±0.0005, re- 
spectively. The line widths and redshifts are consis- 
tent with those measur ed for the CO(J=3— >2) emis- 
sion line (Tacconi et al. I2006D . and are equal for both 
spatial components within the uncertainties. From 
the spatially-integrated emission, we determine a sys- 
temic redshift of z=2.4876±0.0004, which we adopt 
as the nominal value for the system in the following. 
The line parameters for J123707+6214ne and sw cor- 
respond to velocity-integrated emission line strengths 



inates the integrated line emission (Fig. [21 left), extracted 
over a velocity range comparable to the CO(J=1^0) 
map. The maximum signal-to-noise ratio on the sw com- 
ponent is obtained over a narrower velocity range (Fig. 
[21 middle). In this map, the sw component is brighter 
than the ne component, comparable to w hat is seen i n 
the CO(J=3^2) maps of Tacconi et al. (pOOl [2001 . 
Only ~60% of the emission from the ne component are 
seen over this narrower velocity range. We do not detect 
the underlying continuum emission at a 3er upper limit 
of 0.8mJybeam -1 at 1.82mm (rest-frame 520 /mi). 

From Gaussian fitting to the integrated 
CO(J=5^4) line profile (Fig. [2J right), we ob- 
tain S„=4.1±0.8mJy at dv=485±110 kms -1 , cen- 



tered at z=2.4875±0.0006. 



I, 



CO(5-4) 



=2.12±0.51Jy kms" 



This corresponds to 
Fitting the ne and 



of I, 



CO(1-0) : 



=0.180±0.029 and 0.149±0.029 Jy kms' 



i.e., line luminosities of L' co ^ 10 -j=(5.36±0.86) and 

(4.43±0.85)xl0 10 K kms -1 pc 2 . 

We have also detected spatially resolved CO( J— 5— s4) 
line emission at >5cr significance toward both compo- 
nents of J123707+6214 (Fig. [2]). The ne component dom- 

9 Following the nomenclature by Tacconi et al. (2006). 

10 CO(J=l— >0) spectra are Hanning-smoothed. 



sw components individually yields S l „=2.4±0.5 and 
2.1±0.5mJy, dw=467±124 and 432±130 kms -1 , and 
/ C O(5-4)=l-17±0.33 and 0.94±0.29Jy kms -1 , respec- 
tively. We thus derive L^ o j 5 _ 4 ^ = (1.39±0.32) and 

(1.12±0.29)xl0 10 K kms -1 pc 2 , respectively. 

This implies CO J=3^»2/l— >0 line brightness temper- 
ature ratios of r 31 =0.39±0.09 (neJ3 and 0.37±0.10 (sw), 
and CO J=5— >4/l— s>0 line brightness temperature ratios 
of r 51 =0.26±0.07 and 0.25±0.08. The CO(J=3^2) and 
CO( J=5->4) emission lines are clearly subthermally ex- 
cited toward both components (f3i<l and r5i<l). In- 
terestingly, both components appear to have compara- 
ble gas excitation. The ne component is brighter in 
all CO transitions. This suggests that the ne compo- 
nent carries the dominant fraction of the molecular gas 
mass in this system. Based on a ULIRG conversion fac- 

11 We recomputed the CO ( J=3— »2) fluxes based on the data 
presented in Tacconi et al. 120061) by extracting the emission 
over the same (broader) velocity range as the CO(J=l— >0) and 
CO(J=5— ¥&) emission. This yields I, 



CO(3-2) 



=0.63±0.10 and 



0.50±0.10 Jy kms -1 for the ne and sw components, respectively. 
The different methods for extracting fluxes are likely responsible for 
differences between our analysis and an inde penden t CO(J=l— >0) 
study carried out in parallel by Ivison et al. 120111) . 
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Fig. 3. — Channel maps of the CO(J=l— >0) emission in J123707+6214. One channel width is 182 kms -1 (20 MHz; frequencies increase 
with channel number and arc shown at 33.010, 33.030, 33.050, 33.070, and 33.090GHz). Contours are shown in steps of 1<t=57 /iJy beam - 1 , 
starting at ±2<r. The beam size (as indicated in the bottom left) and crosses are the same as in Fig. [T] 



to derive M gas from 



tor a C o=0.8M Q (K kms^pc 2 ) 

-^CO(i-o) (Dowries & Solomon [1998), we determine the 
total molecular gas masses of J123707+6214ne and sw 
to be M gas =4.3 and 3.5xl0 10 M Q [3 i.e., by more than 
a factor of 2 higher than previously found based on the 
CO(J=3— >2) data (scaled to the same aco)> an d cor- 
responding to ~2/3 of th e stellar mass in this system 
(Tacconi et al. 120061 I2008D . 

3.2. Dynamical Structure of the Gas Reservoir 

In Figure [31 maps of the CO( J=l—>0) emission are 
shown in 182 kms -1 wide velocity channels. The emis- 
sion toward J123707+6214ne appears dynamically re- 
solved on — 1.5" (—12 kpc) scales, which may suggest 
that the emission is more spatially extended than in the 
C O(J=3 ->2) line (0.5"±0.2", or 4.1±1.6kpc; Tacconi et 
al. 120061) . J123707+6214sw appears marginally spatially 
resolved in position- velocity space at best, consistent 
with the size measured in CO(J=3— >2) emission within 
the relative uncert ainties (0.9"±0.3", or 7.4±2.5kpc; 
Tacconi et al. [2006). Assuming radii of 6 and 3.7kpc for 
J123707+6214ne and sw, this yields dynamical masses of 
Mdyn sin 2 i=2.9 and 1.5xlO 11 M0 (which we estimate to 
be reliable within a factor of 2). This is about twice as 
high as previous est imate s based on CO(J=3— >2) emis- 
sion (Tacconi et al. 12008ft . and corresponds to gas mass 
fractions of / gas =0.15 and 0.23 for J123707+6214ne and 
sw, respectively. Higher resolution observations are re- 
quired to better constrain how the merger dynamics im- 
pact the CO line profiles and the morphology of the gas 
reservoir, which is necessary to determine more precise 
dynamical masses. 

4. ANALYSIS AND DISCUSSION 

4.1. Origin of the CO Emission 

SMGs are commonly associated with heavily obscured 
starbursts. J123707+6214 is a particularly insightful ex- 
ample of this population, as the ne component remains 
undetected at all wavelengths shortward of 3.6 /im (rest- 
frame 1.0 fim). As shown in Figure HI the peak of the 
CO( J=l— >0) emission of the ne component is clearly as- 
sociated with peaks in the mid-infrared (8.0 fim; rest- 
frame 2.3 fim) and radio co ntinuu m (20 cm; rest-frame 
6 cm; see also Tacconi et al. 12008ft . The ne component 
also slightly dominates the radio emission (—55%), which 
suggests that it contributes the dominant fraction to the 
source's SFR. As also shown in Fig. gj J123707+6214sw 

12 A M ilky- W ay-like aco=3.5 Mq (K kms -1 pc 2 ) -1 (e.g., 
Daddi et al. 2010) would increase M gas by a factor of 4.4. 



consists of multiple components in the optical (606 nm; 
rest-frame 174 nm) that a re separated by a few kpc (see 
also Swinbank et al. l2004[ ). This may correspond to mul- 
tiple star-forming clumps, embedded in a more complex, 
extended molecular gas reservoir, but may also reflect 
the high degree of obscuration in this source. 

4.2. Gas Surface Densities and Star Formation 
Timescales 

Estimating that the CO(J=l— >0) emission in 
J123707+6214ne and sw is distributed over 6 and 
3.7 kpc radius regions, the surface-averaged gas densities 
are E gas =3.8 and 8.1xlO 8 M0kpc -2 , i.e., comparable 
to but somewhat lower than estimates based on the 
(assuming the above size estimate s) more co mpact 
CO(J=3->2) emission (Tacconi et al. I2006L l2008h . This 
would be consistent with some of the CO(J=l— >0) 
emission being in a diffuse, low surface brightness 
component. 

Based on th e SFR of 5OO±25OM yr" 1 determined by 
Tacconi et al. (2008J), we derive a gas depletion timescale 
of T dc p=A:/ gas /SFR -160 Myr for J123707+6214. This is 
consistent with but on the high end of what is found for 
other SMGs (for whic h M ga s are inferred from mid- J CO 
lines; e.g., Greve etal.[2005). 

5. CONCLUSIONS 

We have detected spatially resolved CO(J=l— M3) 
and CO(J=5->4) emission toward the z=2.49 SMG 
J123707+6214. We resolve the emission into two compo- 
nents pre viously dete cted in CO( J=3— >2) emission (Tac- 
coni et al. [2006, 2008), which are likely merging galaxies. 
Both components show similar CO excitation properties, 
with moderate J=3— >2/l— >0 line ratios of r3i— 0.38, and 
relatively low J=5— >4/l— >Q line ratios of rsi — 0.25. The 
implied J=5— >4/3— >2 line ratios of r53 — 0.66 are compa- 
rable to those found in other SMG s that show evidence 
for mergers (e.g., WeiB et al. 120051) . On the other hand, 
the low are comparable to those found in massive gas- 
rich star-formin g gala xies with much lower SFRs (Dan- 
nerbauer et al. [2009; Aravena et al. [20Tot ). This may 
suggest that, in addition to the highly-excited gas as- 
sociated with the starburst, J123707+6214 hosts a sub- 
stantial amount of low-excitation gas. 

The CO(J=l— s-0) emission suggest the presence of 
>2.5x more molecular gas than expected if assuming 
a constant brightness temperature from CO(J=3— >2). 
The optically detected merger component (sw) carries 
~45% of the gas mass in this system, suggesting compa- 
rable amounts of gas in both components, with a slightly 
higher contribution coming from the optically invisible 
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component (ne; *~55%). The radio continuum emission 
consistently indicates a comparable starburst strength 
in both components. Assuming that none is substan- 
tially contaminated by an obscured AGN, and given the 
high SFR of 500±250M Q yr -1 , this provides support- 
ing evidence for a ULIRG-like aco in both components. 
The CO(J=1^0) emission in J123707+6214 likely arises 
from the same gas phase detected in the higher- J lines, 
but the CO(J=l— ^0) emission appears somewhat more 
spatially extended. This yields a revised, ^2x higher es- 
timate for the dynamical mass of the system. Also, this 
finding would be consistent with the presence of some dif- 
fuse, low-excitation gas (which may have a higher aco 
than the highly-excited gas). Such a low-excitation com- 
ponent could be associated with gas that is redistributed 
by mechanical energy input from the starburst, or with 
tidal structure in the ongoing, gas-rich merger in this 
system. 

Our findings highlight the importance of observing 
multiple CO lines including CO(J=l— s-0) to determine 
the total molecular gas mass and gas propert ies in SMGs 
(as already acknowledged by Tacconi et al. 120081 in the 
initial observations of this source). Our results are con- 
sistent with those found for other S MGs o bserved in 
C O(J=1 -»0) emission (Hai nline et al. 12006] Carilli et 
al. 120101: Ivison et al. [2010; Harris et al. 120101 ) . which 
commonly show lower CO line excitation than typically 
found in FIR-luminous quasar host galaxies at compa- 



rable redshifts and with comparable gas m asses (e.g., 
Riechers et al. [2001: [200l Weifi et al. l2007l ). This pro- 
vides supporting evidence that both populations trace 
different evolutionary stages of the same massive galaxy 
population, as would be expected in the ULIRG-quasar 
transition scenario proposed by Sanders et al. (|1988[ ) . 

J123707+6214 is a prototypical example of an SMG 
during an early merger stage, found in the peak epoch of 
galaxy formation. Higher resolution, dynamical mapping 
of CO( J=l— 7-0) emission in this intriguing system (and 
others) is desirable to narrow down aco through dynam- 
ical mass measurements over several resolution elements, 
as possible with the full EVLA in the future. A more 
complete census of CO(J=l— >0) observations of SMGs 
will provide the necessary context to interpret the results 
of such investigations. Such studies provide the most di- 
rect means to constrain the gas fraction, total mass and 
evolutionary state of SMGs, which is necessary to bet- 
ter understand the evolutionary path of massive galaxies 
through their most active phases, and to constrain the 
molecular gas mass density of the universe. 
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